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We study long-time properties of reversible reaction-diffusion systems of type A + B «-» C by 
means of perturbation expansion in powers of 1/t (inverse of time). For the case of equal diffusion 
coefficients we present exact formulas for the asymptotic forms of reactant concentrations and a 
complete, recursive expression for an arbitrary term of the expansions. Taking an appropriate 
limit we show that by studying reversible reactions one can obtain "singular" solutions typical of 
irreversible reactions. 
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I. INTRODUCTION 

Behavior of many physical, biological, and chemical 
systems is determined by evolution of a reaction front 
formed between initially separated reactants jl], g]. The 
simplest theoretical model of this phenomenon consists in 
assuming that initially two species A and B are uniformly 
distributed on opposite sides of an impenetrable barrier. 
The barrier is removed at time t = and the two species 
start to mix and react, forming a dynamic reaction front. 
It is assumed that diffusion is the only transport mech- 
anism and that the diffusion coefficients of each species 
are strictly constant, i.e., independent of spatial location, 
reactant concentrations, etc. The problem to be solved is 
to give a detailed description of spatiotemporal evolution 
of this reaction-diffusion system. 

Simple as it is, this system exhibits many unexpected 
features, especially when the reaction is irreversible (A 
+ B — > C). Using a scaling ansatz, Galfi and Racz 0] 
showed that in this case the width of the reaction front 
grows asymptotically as t a with surprisingly small value 
of the exponent a — 1/6. Renormalization group tech- 
nique , dimensional analysis ^ || , and extensive com- 
puter simulations JtJ were then employed to demonstrate 
that a = 1/6 only above the critical dimension d c = 2 
and for d < d c one must take into account fluctuations 
of reactant concentrations. It was also shown @, |) that 
the case where one of the reactants is immobile (Da = 
or Db = 0) belongs to a separate universality class, with 
a = 1/2. Short-time perturbation flTo| and numerical 
jll| analysis revealed that the reaction front can move in 
a nonmonotonic manner. The quasistatic approximation 
H [l2| was used to find a detailed description of con- 
centrations of species A, B, and C outside the reaction 
zone H [l3|, |l4| and the case of equal diffusion coefficient 
was treated rigorously These theoretical results are 
in full agreement with experiments |l6|, 17, |l8[ [lj| |2^| 
and were generalized to several more complex reactions 
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In reality, however, most of chemical reactions are re- 
versible. In spite of this, reversible A + B ^ C reaction- 
diffusion processes did not attract so much attention. 
This should be probably attributed to the fact that re- 
versible reaction-diffusion processes for a long time were 
not supposed to exhibit any "anomalous" properties, es- 
pecially after Chopard et. al 63 showed that (a) the front 
width of a reversible reaction asymptotically scales with 
time as if the process was governed only by diffusion 
(w(t) oc t 1 / 2 ) and (b) the fluctuations do not modify 
the scaling exponents even in one-dimensional systems. 
However, the problem of giving a detailed description of 
spatiotemporal evolution of reversible reaction-diffusion 
systems remained open. 

This issue was recently considered by Sinder and Pel- 
leg in a series of papers [jl4], [50j. They restricted 
their analysis mainly to systems with very small back- 
ward reaction rate g and found that in this case con- 
centrations of species A, B, and C are practically the 
same as those observed in strictly irreversible reactions 
(g = 0) everywhere except in a very narrow reaction 
zone. They confirmed the result of Ref. ]2£| that there is 
a crossover between intermediate-time "irreversible" and 
long-time "reversible" regimes. They also showed that 
in contrast to irreversible systems, the effective asymp- 
totic reaction rate R can have two maxima and there 
can be even a region where it is negative. Moreover, 
they presented strong arguments in support of a conjec- 
ture that reversible reaction-diffusion processes belong to 
two distinct universality classes. One of them consists of 
systems with immobile reaction product C and asymp- 
totically immobile reaction front, while all other systems 
belong the other universality class. 

The aim of our paper is to work out a new tool for anal- 
ysis of arbitrary reversible reaction-diffusion systems — 
expansion of concentrations as series in 1 jt (where t de- 
notes time) — and to use it to find exact solutions of the 
problem at least for some particular combinations of con- 
trol parameters. This kind of approach was already suc- 
cessfully applied to reaction fronts at very short times 
JlOt - To our knowledge, however, such a technique has 
not been applied to the long-time regime. 

The structure of the paper is as follows. In section II 
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we give th e prec ise definition of the problem to be solved. 
In section [II A we show how to apply the perturbation 
expa nsion to reversible reaction-diffusion systems. In sec- 



tion IIIB 



we employ this technique to study the case 
where diffusion coefficients of all species are the same. 
In particular, we present (a) explicit forms of the long- 
time reaction rate R and the concentrations of species A, 
B, and C; (b) a reccursive formula for any term of the 
expansions; (c) detailed analysis of the cross over to the 
long-time, reversible regime. In section III C we analyze 
two particular cases where some aspects of the system 
behaviour can be studied analytically. Finally, section 



[V is devoted to discussion of results. 



II. MODEL 

We shall consider the following system of nonlinear 
partial differential equations || 

da(x,t) D /\(^) _ R{Xit) (1) 



dt 



dx 2 



db(x,t) d 2 b(x,t) 

= D B — — R(x,t) (2) 



dt 



dx 2 



dc(x,t) Dc d 2 <^) +RM (3) 



dt. 



dx 2 



where 



R(x, t) = ka(x, t)b(x, t) — gc(x, t). 



(4) 



Here a(x,t), b(x,t), and c(x, t) are concentrations of 
species A, B and C, respectively, Da, Db, and Dc are 
their respective diffusion coefficients, R(x, t) denotes the 
effective reaction rate, while k > and g > are the for- 
ward and backward reaction rate constants, respectively. 
The initial conditions to ([!]) - (ji|) read 

a(x,0) = a H(x), b(x, 0) = b H (-x), c(x, 0) = 0, 

(5) 

where ao,&o are initial concentrations of species A and 
B, respectively, and H(x) is the Heaviside step function, 
which is for x < and 1 for x > 0. Note that in 
Refs U, |29[ |3fj | R was called a "refined reaction rate" 
and denoted as R r . Our main goal is to give a detailed 
description of the long-time solutions of ([!]) - (^). 

As was pointed out in Ref. p£|, by measuring length, 
time, and concentration in units of ^/DA/kao, 1/kao, 
and do, respectively, our original problem can be reduced 
to the one with 



D A 



«o 



1. 



1. 



(6) 



We shall adopt these particular values in our further anal- 
ysis, which will leave us with four independent control 
parameters: g, bo, Db, and Dc- 



III. THE ASYMPTOTIC EXPANSION 

A. Formalism 

The analysis performed in Refs. jli], |2^, ^9) shows that 
in the asymptotic, long-time limit solutions to equations 
(0) " © ^ a ' te 011 a scaling form 



a(x,t) w SA(x/Vi), 

b(x,t) w S B (x/Vt), 

c(x,t) w S c {x/Vi), 

R{x,t) w t- l S R (x/Vi). 



(7) 
(8) 
(9) 
(10) 



We therefore assume that in the long-time limit all func- 
tions involved in (|l|) - (Q) can be expanded as series 
in r = 1/t, with coefficients being some functions of 
£ = xj\ft. We thus write a(x,t), b(x,t), and c(x,t) as 



Using 



XVMn(0. (11) 
n=0 
oo 

X>«B„(e), (is) 

oo 

Y,r n C n (0- (13) 
we conclude that R(x,t) can be expressed as 

oc 

R{x,t) = Y j r n 1l n {0, (14) 



a(x, t) 
b(x,t) 
c(x, t) 



where 



K n (0 = -gCnft +J2M0®n-j(0- (15) 

3=0 

Substituting (pj) - Q into (g) - (|) and collecting 
coefficients at t° we find that 



tto(0 = -4>(OBo(0-flCb(0 = 0. 



(16) 



This equation expresses a fundamental property of the 
system: in the long-time limit it tends to a local chemi- 
cal equilibrium (the forward and backward reaction rates 
become asymptotically equal) p8[| . Collecting now coef- 
ficients at r n+1 , n = 0, 1, . . ., we arrive at 



i 

f 2^ 





= 


(17) 


ftn+l(£) 


= 


(18) 




= 


(19) 
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where we used a short-hand notation /'(£) = df/d£, 
/"(£) = d 2 f/d£_ 2 . In the lowest order (rc = 0) we thus 
find 



+ = 
^'(0 + k«o(0 = ^i(0 



^'(0 + 2^(0 = -^(0- 



(20) 
(21) 
(22) 



Equations (|16|), (|2C|) - ( |22j ) constitute a system of four 
equations for four unknown functions: Ao(^), Bq(£), 
Cq(£), and 7?-i(£). These important functions control the 
asymptotic (long-time) properties of the system and can 
be readily identified with the scaling functions employed 
in relations (Q) — (JlC }. However, owing to a nonlinear 
form of the equation (16|), the explicit form of these func- 
tions can be found only in a few particular cases discussed 
below. 



B. The case of equal diffusion constants 

1. Asymptotic solution 

Following (^J) we will now assume Da = Db = Dc = 
1, k = 1 and a = 1. The only free parameters of the 
problem are thus bo and g. Upon adding (||) to ([!]) and 
(^) we arrive at two diffusion equations with well known 
solutions 

(23) 

b(x, t) + c(x, t) = ^b crk(-x/V4:t). (24) 
This immediately leads to 

AM) _ .<a-, + (26) 
m) _ -« ( fl-, + ^it? (26) 

Co(?) = crfc({/2) + Mrfc(-|/2) + 2g - ^VM^ 7) 



a(x,t) + c(x,t) — — erfc(x/v4i), 



5 6 (l+.g + 6 )cxp(-C 2 /2) 



27r[A(C)] 3 / 2 
A n (0 = B n (0 = -C n (0, n > 1 
where we denoted 

1 



(28) 
(29) 



= -[a erfc(e/2)-6 erfc(-e/2)] (30) 



A(0 = [<!>(£) - ff ] 2 + 2 3 erfc(£/2). 



(31) 



Note that 7£i(£) can diverge to infinity. This can hap- 
pen if and only if A(£) — > which, in turn, occurs only if 



<&(£) — > and g — ► 0. Since <&(£) decreases monotonically 
from ao to — &o, equation <&(£) = has a unique solution, 
which will be denoted £/. 

In the limit of an irreversible reaction we find 



lim Ao{£) 

g^0 



lim B (0 

g^rO 



0, 



0. 



limC„(0 = / l & oerfcK/2), f < 

9^0 



limtti(0 

3^0 



6oexp(-^/ 4 ) 



(32) 



(33) 



(34) 



(35) 



V^erfc(C//2) 

where 5 is the Dirac's delta distribution. These relations 
are in full agreement with general formulas derived in 
Ref. jl3| for the case of a strictly irreversible reaction, 
9 = 0, and in Ref. Q for a — b and g — » 0. Both 
species, A and B, become effectively segregated at 
We can thus identify £/ as the position of the reaction 
front. Beyond this point the reaction rate tends to and 
the concentrations a(x,tY b{x,t), and c(x,t) satisfy dif- 
fusion equations (|l|) - (j3|) with R(x, t) = 0. However, 
the existence of a singularity in the above formulas does 
not mean that the reaction should be asymptotically re- 
stricted to a single point! Recalling that £ = x/y/t we 
conclude that this singularity indicates that the reaction 
must be restricted to a region much narrower than it 
also suggests that for g — and xj\ft » we should try 
and take into account higher-order terms of expansions 

(0) - ©■ 

In the oposite limit of an infinitely large backward re- 
action rate g we find 

lim Ao{0 = o erfc (£/ 2 ) ( 36 ) 



lim B (0 



lim C„(0 

g-*oo 



1 



6 erfc(-e/2) 



lim 7£i(£) = 0. 



(37) 
(38) 



These equations express the fact that for large backward 
reaction rates g if a particle C is being created as a result 
of a forward, A + B — ► C reaction, it is being immedi- 
ately converted back into a pair A-B; consequently, the 
concentration of particles C tends to and the concen- 
trations of particles A and B evolve as if there was no 
reaction at all. 



2. Recursive formula 



To find the remaining terms of expansions (11)~(|R 
we employ relations ([l5|), (|lS|), (|22|), and (p9|), arriving 
at a recursive formula 



c„+i(0 - 



(39) 
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where S n {Q = S"=i Cj(£)Cn+i-j(£) and n > 0. An im- 
portant feature of this relation is that it allows to express 
the (n+l)-th term of the expansion directly as a function 
of already determined, lower-order terms. Consequently, 
together with the explicit form of the zcroth order terms 
given in (^BJ) - (|28l ) , equations (|29|) and (^9|) enable one to 
calculate (at least in principle) an arbitrary term of the 
expansions ([n]) - (|1|) analytically. In practice, however, 
the complexity of the appropriate formulas grows rapidly 
and even using computer-algebra systems it is very dif- 
ficult to determine C n (£) for more than a few smallest 
values of n. Notice also that for £ = £f and g — > the 
denominator of (J3S|) goes to 0. 



3. Crossover from "irreversible" to "reversible" reaction 
fronts 

Let t* denote the time when the system enters the 
asymptotic, long-time time re gim e. For times t t* 
each sum in expansions (^) - ( p"3| ) should be dominated 
by its lowest-order nonvanishing term, while for t%t* the 
system should behave as if the reaction was strictly irre- 
versible (<7 = 0) JlJ; H^]. This cross-over time can be esti- 
mated from a relation -Ao(£/) ~ ^4i(£/)A* ■ Using (22), 
©, and © we find A(£) = -C x (0 = fti (0/^/5(0. 
In the limit 5 — > there is thus .Ao(£/) K ana - 
.4i(£/) oc Therefore in this limit we have 



t* oc 



-3/2 



(40) 



Note that Chopard et.al. Q proposed a different re- 
lation, t* oc However, their conjecture was based 
on numerical analysis of a reaction-diffusion system with 
a very specific choice of system parameters: clq = bo, 
Da = Db-i and Dq = 0. In other words, they inves- 
tigated only symmetric systems with immobile particles 
C. Moreover, they assumed that the width of the reac- 
tion front can be identified width the width of the profile 
of particles C, wc{t) = J x 2 c(x,t)/ J c(x,t)dx. Actually 
this is acceptably only when both particles C and the re- 
action front center are immobile (Dc — 0, Xf(t) ~ 0), a 
condition implicitly satisfied in their simulations. For 
Dq 7^ or Xt(t) 7^ we expect that asymptotically 
w c (t) oc t 1 / 2 , while w(t) oc i 1/6 (at least for g = [§), 
so w c (t) cannot be identified with w(t). From this point 
of view the case studied numerically in Ref . [^8| belongs 
to a separate universality class. To the same conclu- 
sion, though on different grounds, came Sinder and Pel- 
leg ps| , who investigated systems with immobile reac- 
tion product (Dq — 0). They found that asymptotically 
w(t, g) oc g 1 / 2 ^/ 2 for systems with the moving reaction 
front (x f ^ 0) and w(t,g) oc g^H 1 / 2 if x f ~ 0. There- 
fore, on taking into account that for small g and inter- 
mediate times t one expects w(t,g) oc g°t 1 / 6 , we immc- 
diatelly arrive at the conclusion that for asymptotically 
immobile reaction fronts studied in Ref. J2§] there should 
be t* oc while for systems with Dc = 0, g <C 1, and 



a mobile reaction front (xt 7^ 0) the crossover time is 
given by ©. 

We verified the validity of relation ( |40| ) for systems 
studied in this section by computer-assisted analysis of 
An(£f) for the fully symmetric case ao = &o = 1. It 
indicates that 



A*(£/) oc (- 



-l)»g«- 



3n)/2 



(41) 



for all n < 5. Most probably this relation continues to be 
true also for higher-order terms. This would mean that 
for t ^> t* , as expected, expansions (|ll| ) - (|l3| ) will be 
dominated by their first terms, while for fct* they are 
divergent at £/. 

Using ( f4(i| ) we can estimate the width w and 
height h of the reaction front at the cross-over time. 
The former is defined as a square root of J_ (x — 
Xf) 2 R(x,t)dt/ R(x,t)dt, while h = R(xf,t) (here 
Xf w VtCf is the exact location of the reaction front, 
see Ref. j| |§). Usin S © and © we nnd 



^*) ex (t*rvv?cx(o- 2 / 3 . 



(42) 
(43) 



These are the scaling relations derived by Galfi and Racz 
for the case of a strictly irreversible reaction, g = and 



C. Other cases with rigorous solutions 

1. The limit |£| — ► oo, or the tails of the distributions 

Let us now consider the general case of arbitrary values 
of parameters 6 , Dg, Dq and g in the limit of t — > 00 and 
|a;| — > 00 such that l^l/v 7 ^ = |£| — * 00 ■ F° r sufficiently 
large £ we may expect that the concentration of particles 
B will be very close to its original value bo- Substituting 
B {£) ~ b in @ - © we find that K^) « and 

Cb(0 « feoff-Mo (0 « rj + erfoU/J^] (44) 



where 



+ _ 9Da + &o-Dc 



9 + bo 



(45) 



and rj + is an integration constant independent of £. Note 
that, as might be expected on physical grounds, D~^ s lies 
between Da and Dc- 

Similarly, in the opposite limit £ — * —00, we assume 
Aq (£) ~ ao = 1 and find 



C (0 « ff-'BoCO « rf erfc -^NAD~ S 



where 



d; 



gD B + Dc 
5+1 



(46) 



(47) 
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and T]~ does not depend on £. We thus see that suffi- 
ciently far away from the reaction region all the species 
diffuse with an effective diffusion coefficient (for 
£ — > oo) or (for £ — > — oo). Note that D^, 

D~ s — > as g — > 0, in agreement with the findings 
ofRef. Q. 



Reaction front at £ = /or one or too vanishing 
diffusion constants 



Equations ( |2l| ) - (22) immediately imply that if or 
vanishes then 



fti(0) = 0. 



(48) 



Thus either has a local minimum at £ = or it 

attains negative values in the vicinity of £ = 0. Both pos- 
sibilities have actually been observed in numerical simu- 
lations carried out in Ref. E9I. 



IV. CONCLUSIONS 

We have applied perturbation analysis to reversible 
reaction-diffusion systems with arbitrary values of control 
parameters. Using this technique we obtained a system 
of four equations completely governing the long-time be- 
haviour of the system. We then concentrated on several 
cases where analytical results can be derived. 

In particular, we found the complete solution of the 
problem for the case where the diffusion coefficients of 
all species are the same. Its most interesting feature 
is the limit of vanishingly small backward reaction rate 
constant g. In this limit, as expected, the solutions be- 
come singular at a point which can be identified with 
the reaction-front center. Our method enables one to 
analyze explicitly how these singularities, especially the 
Dirac's delta function in the expression for the local re- 
action rate 72-1, appear in mathematical formulas. This 



holds out hope that it will be possible to construct a uni- 
fied theory of reversible and irreversible reaction-diffusion 
systems. It should be also noticed that our asymptotic 
solution is also valid for g — even though in this case 
all higher-order terms of the expansion diverge. 

We also found a general recursive formula for any term 
of the expansion as a function of all lower-order terms. 
Unfortunately, complexity of expressions thus obtained 
grows very rapidly. Nevertheless we believe that this 
formula opens a new way of analyzing the evolution of 
reaction-diffusion systems at arbitrary times. 

Using the information about the first correction to 
the asymptotic solution we showed that the anomalous 
properties of irreversible reaction-diffusion systems can 
be studied by taking a suitable limit in the formulas 
obtained for reversible systems. A potential advantage 
of this approach is that mathematical description of re- 
versible reaction-diffusion systems is more regular, and 
hence more amenable to rigorous analysis. 

We also studied the crossover time t* between 
intermediate-time ( "irreversible reaction" ) and long-time 
("reversible reaction") regimes. We proved that in the 
case of equal diffusion coefficients t* scales with g as 



-3/2 



This conclusion agrees with results obtained by 



Sinder and Pelleg |29| and disagrees with those obtained 
by Cornell and Droz This discrepancy can be easily 
understood if one notices that Cornell and Droz studied 
only systems where the width w c of the concentration 
of species C grows in time as t 1 / 6 , while our study was 
performed for a system where w c oc i 1 / 2 . 

The main message coming from our work may be sum- 
marized as follows: (a) reversible reaction-diffusion sys- 
tems are more amenable to rigorous treatment than ir- 
reversible ones; (b) it is possible to investigate the more 
difficult irreversible reaction-diffusion systems by taking 
an appropriate limit in the reversible ones; (c) investiga- 
tion of reversible reaction-diffusion systems is interesting 
not only per se, but constitutes an alternative technique 
of analyzing irreversible systems. 
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